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Within the mode-coupling theory for ideal glass-transitions, the mean-squared displacement and the correla- 
tion function for density fluctuations are evaluated for a colloidal liquid of particles interacting with a square- 
well potential for states near the crossing of the line for transitions to a gel with the line for transitions to a 
glass. It is demonstrated how the dynamics is ruled by the interplay of the mechanisms of arrest due to hard- 
core repulsion and due to attraction-induced bond formation as well as by a nearby higher-order glass-transition 
singularity. Application of the universal relaxation laws for the slow dynamics near glass-transition singularities 
explains the qualitative features of the calculated time dependence of the mean-squared displacement, which are 
in accord with the findings obtained in molecular-dynamics simulation studies by Zaccarelli et. al [Phys. Rev. 
E 66, 041402 (2002)]. Correlation functions found by photon-correlation spectroscopy in a micellar system by 
Mallamace et. al [Phys. Rev. Lett. 84, 5431 2000)] can be interpreted qualitatively as a crossover from gel to 
glass dynamics. 
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I. INTRODUCTION 

The mode-coupling theory for ideal glass transitions 
(MCT) is based on closed equations of motion for the corre- 
lation functions of the density fluctuations p^ of wave vector 
q, $ q (t) = (pi(*)P?)/<|P?l>. q = \q\ QH- The static structure 
factor S q enters these equations as input; it is assumed to be 
a smooth function of the control parameters like density p or 
temperature T. The equations of motion exhibit bifurcations 
for the long-time limit of the correlators, f q = Um,^ c(j), ? (r), 
which are referred to as glass-transition singularities. Only 
bifurcations of the cuspoid family can occur in the MCT equa- 
tions Qyll, i.e., singularities of the class A/, / ^ 2, which are 
equivalent to the bifurcations in the real roots of real polyno- 
mials of order I [4]. The generic singularity when changing a 
single control parameter is the A 2 also called fold. In the most 
important situations, it deals with the transition from a liquid, 
characterized by f q = 0, to an idealized glass, characterized 
by f q > 0. The quantity f q is the Debye- Waller factor for the 
arrested amorphous structure. For parameters near a glass- 
transition singularity, slow dynamics emerges with subtle de- 
pendence on time and control parameters. This dynamics is 
proposed by MCT as the explanation for the structural relax- 
ation in glass-forming liquids. The universal laws for this dy- 
namics can be obtained by asymptotic expansion of the equa- 
tions of motion as was demonstrated comprehensively for the 
hard-sphere system (HSS) |5j|6J]. The glass transition for the 
HSS has been studied experimentally by dynamic hght scat- 
tering for sterically stabilized hard-sphere colloids 00,0]. 
The successful analysis of the data within the MCT frame pro- 
vides strong support for the theory 1 10]. 

It is known from studies of so-called schematic models, 
that there may emerge also higher-order singularities from 
MCT like A3 and A4 Hill . The most significant feature of 
the dynamics near an A/ with I ^ 3 are logarithmic decay 
laws, where detailed properties have also been worked out 
in full generality 11211 . There is a variety of data indicating 
that logarithmic decay laws occur in some glass-forming liq- 
uids Gj 0O E E3 Generically, one has to vary 



two or three control parameters, respectively, in order to ap- 
proach these higher-order singularities. It was discovered only 
recently that the MCT equations for simple systems imply 
the existence of an A3-singularity if a hard-sphere repulsion 
is complemented by a short-ranged attraction shell lfl9l Eol . 
The A3 is the endpoint of a line of A2-singularities describing 
glass-to-glass transitions in the parameter plane spanned by 
the packing fraction cp and the effective attraction strength Y. 
At this line there occurs a transition from a glass caused by 
the cage effect due to the strong repulsion to a glass caused by 
bond formation due to the dominant role played by the attrac- 
tion. This transition line extends to low packing fraction and it 
was argued to be related to the gel transition there 1 20| . There- 
fore, this line shall be referred to as gel line in the following 
for the sake of brevity. There is a second transition line that 
extends to the known transition of the HSS if Y tends to zero. 
For brevity, this line shall be referred to as glass line in the fol- 
lowing. The glass line terminates transversally at the gel line 
forming a line crossing in the glass-transition diagram. The 
liquid dynamics close to this crossing shall be studied in this 
paper. 

The existence of a crossing point depends on the attraction 
to be sufficiently short-ranged. If the range 8 of the attrac- 
tive potential increases above a critical value, the glass-glass 
transition line and the A3-singularity vanish. This happens 
in an A4-singularity as was demonstrated first for the sim- 
ple system of particles interacting via a square-well poten- 
tial 12111 . The topological singularities A/ are robust against 
parameter variation. It was shown explicitly for a variety 
of cases that various interaction potentials or approximation 
schemes for the static structure factor yield the same quali- 
tative results In this paper, the square-well 
system (SWS) shall be used as model for the quantitative 
work. Systems with short-ranged attraction can be realized in 
colloid-polymer mixtures, where the polymer induces a deple- 
tion attraction 112 511 . Such systems are well under control ex- 
perimentally and have established thermodynamic phase be- 
havior [26]. Logically disconnected from the appearance of 
higher-order singularities, MCT predicts a subtle reentry phe- 
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nomenon for the glass transitions in such systems fl^l which 
can be related to the variation of the static structure factor 1 2 1 ] . 
Starting in the glassy state of the HSS and increasing the at- 
traction, the glass is melted for sufficiently small range of the 
attraction. Upon further increasing the attraction, the system 
arrests again. This reentry phenomenon is now firmly estab- 
lished by experiments in colloidal systems J27l I2I and by 
molecular dynamics simulation I15ll28ll29ll30ll . 

The scenario suggested by MCT for the A2-singularity has 
been applied successfully to analyze experiments and results 
of computer simulations |31]. It was also applied to systems 
where both glass and gel transitions occur Il4 Il5l l3(ll . For 
the dynamics near higher-order singularities, detailed predic- 
tions for logarithmic decay and subdiffusive power law in the 
mean-squared displacement (MSD) have been worked out for 
the SWS 13211 . Indications of logarithmic decay were reported 
fl4ll which are compatible with MCT predictions, and a re- 
cent study identifies both logarithmic decay in the correlation 
functions and a subdiffusive power law in the MSD which 
is consistent with MCT II 611. It is the main objective of the 
present paper to discuss scenarios in the SWS near a cross- 
ing point where the dynamics is influenced by different A2- 
singularities and higher-order singularities at once. There are 
signs of crossing phenomena connected to higher-order sin- 
gularities in recent experiments with photon correlation spec- 
troscopy in a micellar system fl3l fl7l Hill , a suspension of 
PMMA colloidal particles 1281 13311 . a systems of microgel col- 
loids 0H3]> and computer simulation studies Il5ll28ll . 

The paper proceeds as follows. Section UTI introduces the 
equations of motion of MCT. A comparison of the theoreti- 
cal glass-transition diagram with the simulation of Ref. fiHl 
in Sec.[ni]motivates the asymptotic analysis which is outlined 
in Sec. II Vi and applied to the MSD in Sec.|V]and to the corre- 
lation function in Sec. lVll Section fvTTl presents a conclusion. 
The Appendix addresses specific questions arising in the nu- 
merical determination of the glass-transition singularities. 



II. EQUATIONS OF MOTION 

All equations of MCT are based on the equations of 
motion for the normalized density correlators § q (t) = 
(Pjj( f )Pq) / (\Pq\ 2 ) f° r wave-vector q and its modulus q = \q\. 
When Brownian dynamics for the motion in colloids is as- 
sumed, these equations read flll2ll^. l35ll36ll . 

z q d t ty q (t)+$ q (t) + J m q (t-t')d t ,ty q (t')dt' = 0. (la) 

Here, i q = S q /(Doq 2 ), with Do denoting the short-time diffu- 
sion coefficient. S q = (|ps| 2 ) is the static structure factor of 
the system. The initial condition is § q (0) = 1. The kernel is 
a bilinear functional of the correlators, m q (t) = f q [V,0t(f)], 
with 



r ~, 1 / cPk 
?c,[f] = ■= I 7^ v alfkf\a-k\ 



2 J (2tc) 3 3 < kJ 



(lb) 



and the vertex V specified by 

v q ,i = s i s k s \ q -k\P 4- kc k + q-{q-k)c^_^ jq . (lc) 

The direct correlation function c q is connected with S q by the 
Ornstein-Zernike relation, S q = 1/[1 — pcJ I37ll . 

The long-time limit of the correlation function, f q = 
lim r ^oo(j)( / (f ), can be calculated from an algebraic equation, 
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that displays glass-transition singularities when control pa- 
rameters are varied 0] • 

For the dynamics of the tagged particle density, p q (t) — 
exp[/<7? s (f)], one obtains similar equations for the correlation 
function %(t) = (pi*(t)pi) QS, 



^♦£(0+44(0 + £m q (t-t')M q (t')dt' = 0. 



(3a) 



Here r s (t) denotes the tagged particle position, i s q = 1 /(D s q 2 ) 
with the short-time diffusion coefficient for a single particle, 
denoted by D s , We set D = Do in the following. The kernel 
m q (t) = f q [§(t),§ s (t)] is given by the mode-coupling func- 
tional for the tagged particle motion, 



f:\fj s ] = 



p ^2,^7,2 



(231)3 q- 



q-k\ 



(3b) 



For a tagged particle of the same sort as the constituents of the 



host fluid we can set cf = c q . 



The mean-squared displacement (MSD) of a tagged parti- 
cle, 5r 2 (f) = (\r s {t)-r s (0)\ 2 ), obeys 0, 

5r 2 (t)+D s f m^{t-t')hr 2 {t')dt' = 6D s t . (4a) 



The functional m (0 )(f) = lim^o«^(0 = FMSD[§{t),$ s (t)] for 
the MSD reads 



fMSDlf,?} = 



dk 



(67I 2 ) 



2 . P s k (4) 2 f k fi 



(4b) 



The inverse of this functional determines a characteristic lo- 
calization length r s by r 2 = 1 / TMSD[f,f s }- The long-time dif- 
fusion coefficient D s can be defined by lim r ^oo 8r 2 (t)/t = 6D S 
and yields 0] 



D v 



1 



/■ 00 

D s m^\t)dt. 
Jo 



(5) 



For the equations above, the static structure factor S q is 
required as input, which can be calculated from the interac- 
tion potential after some closure relation is invoked l37ll . For 
the square-well system (SWS), we use an approximate ana- 
lytical solution of the mean-spherical approximation (MSA) 
and a numerical solution to the Percus-Yevick approximation 
(PYA) |2ll]. The SWS consists of N particles in a volume V 
at density p = N/V with hard-core diameter d and an attrac- 
tive well of depth mo and width A. We describe the SWS by 
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three dimensionless control parameters, the packing fraction 
(p = d 3 pn/6, the attraction strength F = uo/{k^T) and the rel- 
ative well width 8 = A/d. The unit of length is chosen to be 
d = I. The unit of time is chosen so that Do = 1 /160. Wave- 
vector space shall be discretized to M grid points with spacing 
Aq = 0.4 and a cutoff q ma * large enough to assert convergence 
of the integral in Eq. dlbl for the long-time limit. The proce- 
dures for the numerical solution of Eqs. Q to (|3J have been 
outlined previously 1 32, 38, 39]. Asymptotic laws close to the 
singularities are presented in the appendix that allow for ac- 
curate and fast determination of both A3-endpoints and glass- 
glass transition points. 



III. GLASS-TRANSITION DIAGRAMS 
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FIG. 1: Glass-transition diagram for the SWS using the structure 
factor within MSA. Five cuts through the three-dimensional diagram 
are shown for constant well widths 8 as curves for attraction strength 
r versus packing fraction (p. All curves start at the limit of the HSS 
for r = as indicated by the arrow. For 8 = 0. 1 17 and 0.06 the curves 
(p c (r) vary smoothly as T is increased. The line 8 = 8*= 0.04381 
hits the A4-singularity (*). Curves for 8 < 8* exhibit a crossing point 
(o) and an A3-endpoint singularity (O) as demonstrated for 8 = 0.03 
and 8 = 0.02, where part of the glass-transition line has been erased 
to avoid cluttering the figure. 

The three-dimensional control-parameter space for the 
SWS can be examined by considering cuts through the set 
of glass-transition singularities for constant 8. In each plane 
the transition points are calculated by finding the bifurcation 
points of Eq. (0. Figuref^displays the singularities for several 
cuts. The glass-transition diagram is organized around the A4- 
singularity (*) at ((p*,r*,8*) MSA = (0.5277,4.476,0.04381). 
From there emerge for 8 < 8* both the line of A3-endpoints, 
((p°(S),r°(8)), and the line crossings, ((p°(8),r°(8)), sepa- 
rating glass transitions for r < F° from gel transitions for 
r ^ r°. The line of gel transitions extends beyond the crossing 
point into the arrested state as glass-glass-transition line and 
terminates at the A3-singularity. For 8 > 8*, glass- and gel- 
transition lines join smoothly as seen for 8 = 0.06 and 0.1 17. 



ForS< Sreentry the lines of glass transitions display the reentry 
phenomenon discussed above. At 8 = 8 ree ntiy this reentry dis- 
appears I24H . When using the analytical result for S g in MSA 
we get 

Sreenty = 0- 1 17 ' while for the PYA 

one finds the larger 

value 8^^ = 0.145. To assure that the smaller value for the 
MSA is not caused by the expansion in 8 used for the cal- 
culation of S q , we determine 8 ree ntry again, this time solving 
the MSA numerically. This yields 8^^ = 0.112. There- 
fore the deviation between the MSA and PYA result has to be 
understood as a difference in the way the closure relations in- 
corporate the subtle changes in S q that lead to the reentry as 
explained earlier I21B . 



r 




FIG. 2: Glass-transition diagram for the SWS at 8 = 0.03 (full lines) 
together with isodiffusivity lines for Dq/D s = 10 5 , 10 7 , 10 10 (dashed 
line, from left to right) based on the structure factor using MSA. 
The A3 -singularity is indicated by a circle (O) and a crossing point 
by a diamond (o). On the isodiffusivity lines, states are marked for 
T = 1.67 (+), 5.50 (•), and 6.33 (■). The dotted lines with the 
shaded circle as endpoint show the glass-transition singularities for 
8 = 0.03 based on the structure factor using PYA rescaled in T by a 
factor 5.88 to match the crossing point. 

For the discussion of the crossing we choose the cut 8 = 
0.03 from Fig. \l\ which is shown in Fig. [2] as full line. The 
ratio of the diffusivity D s compared to the short-time diffu- 
sion coefficient D can be used to characterize the distance of 
a chosen state to the liquid-glass-transition line. The dashed 
lines in Fig. [2] show states for constant Dq/D s with D s de- 
fined in Eq. 10. These lines are plotted for the cut 8 = 0.03 
also using the MSA for the evaluation of the structure fac- 
tor. These isodiffusivity lines can be interpreted as approxi- 
mations of the liquid-glass-transition line. They also display 
the reentry phenomenon as discussed above. The liquid-glass- 
transition line follows closely the isodiffusivity curves but is 
separated further from them around the crossing point. This 
indicates the influence of more than one singularity on the 
dynamics in that region. If the PYA instead of the MSA is 
used to calculate the structure factor input, the dotted lines of 
liquid-glass- and glass-glass-transition curves are found. The 
result for both closure relations can be matched reasonably 
at the crossing point by only rescaling T by a factor of 5.88. 
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The agreement for the almost horizontal gel-transition lines 
is less satisfactory but the glass-transition lines almost fall 
on top of each other. As noted in the preceding paragraph, 
the reentry is more pronounced for the result using the PYA 
than for the MSA. The different packing fractions at the cross- 
ing are (p^SA = 0.5364 and cpp YA = 0.5362, while the differ- 
ence in the location of the A3 -singularities is slightly larger, 
<Pmsa = 0.5449 and cp PYA = 0.5456. 
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FIG. 3: Endpoints (O) and crossing points (o) for the SWS in PYA 
for 8 = 8*, 0.04, 0.035, 0.03, 0.025. The crossing points based on the 
MSA can be scaled on top of the PYA result by a 8-dependent prefac- 
tor, T PYA = y(8) T MSA with y(8) ~ 0. 1 + 2.34 8. Crossing points and 
endpoints based on the MSA are shown by filled symbols. The inset 
shows the difference in (p between crossing points and endpoints for 
increasing 8* — 8. Results for the PYA and the MSA are shown by 
open and filled symbols, respectively. The dashed curve displays the 
fit cp* — q>° = 45(8* -8) 2 . 
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FIG. 4: Results for the SWS for 8 = 0.03. Triangles (A) mark the 
isodiffusivity curves from the simulation in 1 15] from left to right 
for D S Q /D S = 2 • 10 2 , 2 • 10 3 , 2 • 10 4 , 2 • 10 5 , respectively. Open 
triangles A indicate the extrapolation of the diffusivity data 1191 . 
Crosses ( x ) show the isodiffusivity curve for DJJD S =2.4- 10 2 from 
the simulation of the monodisperse system [ 29] . Dotted lines are 
guides to the eye for the data from MD simulation. Dashed lines 
indicate the data for melting, freezing and solid-solid binodal to- 
gether with the solid-solid triple point (■) and critical point (•) 
from 14(31 . Full lines are theoretical calculations using the PYA 
structure factor for liquid-glass transitions, the glass-glass transition 
with endpoint A3 (O) and the respective isodiffusivity curves for 
D S Q /D S = 2 • 10 2 , 2 • 10 3 , 2 • 10 4 , 2 • 10 5 (from left to right). The arrow 
labeled HSS indicates the limit of the hard-sphere system from 1 15]. 
The MCT results are based on the PYA and the control parameters 
(p PYA and T PYA are transformed by (p MD = 2.25 (p PYA - 0.5747 and 
r = 2.85 r PYA to match the isodiffusivity curves from the simu- 
lation. 



Figure|3]shows the A3 -singularities and the crossing points 
when using S q in PYA (empty symbols). Matching the cross- 
ing points from the result using the MSA, cf. Fig. [2 again 
by multiplications in Y, yields good agreement in (p° for all 
values of 8. After the transformation, the A3 -singularities for 
a given 8 differ in T by 5% and less, while the deviations in (p 
are comparable to those found for the crossings. It should be 
noted that all endpoints are found at roughly the same attrac- 
tion strength, T w 0.9, whereas the crossing points move to 
higher F as the well width is decreased. At the A/t-singularity, 
the endpoint absorbs the crossing point, and the difference 
cp* — cp° approaches zero in a minimum. Therefore, cross- 
ing point and endpoint separate from each other quadratically 
when close to the A4-singularity. This is demonstrated in the 
inset of Fig. [5] for the results using both MSA and PYA as 
input, respectively. 

One cannot expect a theory for a singularity to predict ac- 
curate numbers for the control parameters of the singulari- 
ties. For that reason the distance from the singularity should 
be used for a comparison of the theoretical results with data 
from experiments or computer simulation. The isodiffusiv- 
ity curves in Fig.[2]motivate a comparison between MCT and 
computer simulation based on the ratio Dq/D s [15]. Figure. |4] 
shows that an acceptable fit of data for the diffusivity in 11511 



and the theoretical data calculated using the structure factor 
evaluated in PYA is achieved by keeping the well width fixed 
at 8 = 0.03 and scaling the axis of the inverse temperature 
by r MD = 2.85 r PYA . This preserves the limiting case of the 
HSS as done above for the comparison of PYA and MSA, cf. 
Fig. [2] Trying to match reasonably at least the two curves 
with the highest ratio of Dq/D s , the packing fraction has to be 
taken (p MD = 2.25 cp PYA - 0.5747 in order to keep a value for 
HSS of (p£j SS = 0.586. This is consistent with the diffusivity 
data and experiments done in colloids 1 8, 9]. The pref actor of 
2.25 seems somewhat large and it is already seen in Fig.|4]that 
this overestimates the differences in <p further from the singu- 
larities. But taking the diffusivity data for granted, this large 
prefactor is required. A modification of the third coupling pa- 
rameter 8 was not necessary in the fit. 

Figure ^demonstrates a reasonable fit between theory and 
data starting from the HSS and extending up to the crossing 
point. For the gel-transitions, there are not enough data avail- 
able to make a definite statement. For this high values of T 
it is also difficult to obtain accurate values for D s with good 
statistics from the simulation 1 15]. These points are only fit- 
ted qualitatively in Fig. |4| An extrapolation of the diffusiv- 
ity data was used in Ref. 1 16] to determine the open triangles 



5 



that represent a different estimate for the liquid-glass transi- 
tion line. These points agree well with the transformed the- 
oretical curves but tend to deviate closer to the crossing. A 
comparison of the fit in Fig. |4] which uses the PYA for the 
theoretical curves, with Fig. |2] indicates that using MSA for 
the structure factor would also properly fit the data from the 
HSS limit up to the crossing but would be worse than PYA 
for the gel line. The indication of the A3-singularity in Fig. [4] 
has to be understood as an extrapolation of the transformation 
scheme outlined above. A slight reservation has to be made 
since the simulation data refer to a binary mixture while the 
present theory deals with a monodisperse system. However, 
comparing the data from the simulation of the monodisperse 
case 1 29] indicated by crosses in Fig. |4] with the ones for the 
mixture, the isodiffusivity for D s /D s = 2.4 • 10 2 seems to fit 
nicely into the picture. Data for lower Di/D s from 12911 have 
the same trend in Y but apparently do not occur at control pa- 
rameter values for the same diffusivity as extrapolated from 
the mixture. The MD studies were performed using Newto- 
nian dynamics where an appropriate definition of Dq is im- 
possible; the value dy/k-gj /m is taken instead of D Q ' as ref- 
erence which introduces a reasonable microscopic time scale 
Jl5l l29ll . This problem in the definition of the analog of D S Q 
introduces less deviations for larger ratios of the diffusivity 
Dq/D s since only the order of magnitude is important for the 
definition of the isodiffusivity curves. A deviation in \ogD 
would stay the same for both large and small differences in 
logDg — logZ> ( and the result can be more accurate the larger 
the ratio D /D s is. Therefore, putting emphasis on the data 
with high ratios of D /D s is justified. 



The fit in Fig. |4] shows that in general MCT overestimates 
the trend to freezing when coupling parameters are increased. 
This was already found for the HSS 1 8] and a binary Lennard- 
Jones mixture 14 ill . Yet, for a Lennard-Jones potential the 
mechanism of arrest is still dominated by repulsion, so the 
control parameter is effectively only density also in that sys- 
tem. For the SWS near the line crossing, necessarily both 
mechanisms of arrest have to be of the same importance and 
the approximation inherent to MCT has to preserve the rela- 
tive importance of both mechanisms. In the case of the SWS, 
MCT has apparently the same tendency in the error for the 
treatment of couplings in cp and Y. The mapping of the theo- 
retical results to higher experimental values of both packing 
fraction and attraction strength is also in agreement with a re- 
cent experimental analysis of a colloid-polymer mixture with 
the theoretical results for the Asakura-Oosawa potential 1 42 ] . 
For the latter work, a qualitatively similar mapping could be 
suggested to match experiments and theoretical predictions. 
By comparison with the data for the phase transitions 1 40] in 
Fig. [4j we recognize that the crossing of lines and the A3- 
singularity are located in the metastable region with respect to 
the solid-solid binodal. The A3 -singularity differs by 4% in (p 
and by a factor of 4.5 in Y from the solid-solid critical point. 



IV. ASYMPTOTIC EXPANSIONS 

For the description of the dynamics at the crossing, asymp- 
totic expansions at the two different singularities shall be 
applied with the separation from the respective singularity 
as small parameter. The separations from an A2- or A3- 
singularity shall be denoted by a and e, respectively. The ex- 
pansions for A2-singularities which are valid for glass-, gel- 
and glass-glass-transition points are taken from Refs. (slBl. 
the expansions for the A3 -singularity are found in I12L 13211 . 
Only those formulas which are needed below are compiled in 
the following. For both A2- and A3 -singularities the expansion 
for the density correlation function can be stated in the general 
form 

<M0 = f q +f q + h q {G{t) + [H(t)+K q G(t) 2 }}, (6) 

where the plateau correction f q and the terms in square brack- 
ets are of next-to-leading order. Neglecting these terms leaves 
the leading order result, tyq(t) = f q + h q G(t), which com- 
prises the factorization theorem of MCT |2], stating that the 
deviation of § q (t) from the plateau f q factorizes into time- 
dependent function G(t) and a critical amplitude h q . This fac- 
torization is violated in next-to-leading order by f q and the 
term K q G(t) 2 with the correction amplitude K q . While the 
general formulas for h q and K q are the same for the ex- 
pansions at both singularities, G(t),H(t), and f q are specific 
for the particular expansion. At an A2-singularity the leading- 
order result is given by the p-correlation function yfl, 

G(?) = v1°Ut('Aa), fa = ro/|a| 1/2a , <^0, (7) 

where the lower signs refer to the weak coupling side of the 
transition. The overall time scale ?o is used as fit parameter. 
For a = 0, the above formula simplifies to a power law as does 
the correction, 

G(t) = (t /t) a , H(t)=K(a)(t /t) 2 ", (8) 

with a function k(jic) 

K(x) = gr(i-3x)-cr(i-^) 3 ]/[r(i-^)r(i-2x)-A,r(i-3x)]. 

(9) 

Here, Y(x) denotes the Gamma function and X is the exponent 
parameter, X = Y(l — a) 2 /r(l — 2a). For an A2-singularity, 
0.5 ^ X < 1, while X — 1 specifies an A3 -singularity. Formulas 
for the parameters Z, and ^ are found in Ref. t aL 
For the MSD, the analog of Eq. © reads KABl 

5r 2 (f )/6 = r c2 ~r 2 - h MSD {G(t) + [H(t) + K MSD G(t) 2 }} , 

(10) 

where only the plateau correction r 2 is again specific to the 
expansion considered. Inserting Eq. (jSJl into Eq. dlOt yields 
the following form for the description of the MSD at the A2- 
transition point |6], 

8r 2 (f)/6 = rf - /imsd (foAHl + [^msd + K(a)](t /t) a } , 

(11) 
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The increase of the MSD above the plateau r£ 2 is given by the 
von Schweidler law, 

5r 2 (f )/6 = rf 2 + / 1MSD (f 1 - [JSTmsd + k(-6)] (f , 

(12) 

with r(l + b) 2 /F(l + 2b) = X. The time scale t' a obeys another 
power-law scaling, t' a = t Q /(B l / b \a\^) ,J = 1 /(2a) + 1 /(2b), 
where the number B is tabulated in Ref. [43]. 

The leading order result for an A3 -singularity is given by 

G(t) = -B\a(t/x) , B= 0-661/tc 2 ], (13) 

where the time scale x is used to match the asymptotic de- 
scription with the solution. The corrections in Eqs. (|6j and 
(II Oi are completed by specifying 

4 

ff(t) = J>ln'(t/T). (14) 

1=1 

The definitions for / ? , r 2 at the A3-singularity and the param- 
eters B, Bj and £1 are found in [ 12]. The solution for the MSD 
at an A3-singularity can be represented in an alternative form 
as a power law 1 32], 

8r 2 (t)/6 = rf(t/zY, (15a) 

with exponent 

x = h MSD B/r c s 2 . (15b) 

The next-to-leading order result implies a correction to the 
exponent 

x' =h MSD (B-B l )/rf. (16a) 

and reads 

5r 2 ( f )/6 = (t/xY{rf - r 2 + b 2 rf\n(t/z) 2 

+ «3 ln(f/x) 3 +a 4 ln(f/x) 4 }. 

Here b 2 = (2r c2 a 2 - a 2 )/(2r c s 4 ), a\ = h MSD (B - Bi), a 2 = 
-hMSD(B 2 +K MSD B 2 ), 03 = -/j M sd#3, and a 4 = -/!msd#4- 

V. RESULTS FOR THE MEAN-SQUARED 
DISPLACEMENT 

Three paths are marked in Fig. [2] for the discussion of the 
dynamics. The first path for T = 1.67 is relatively far from 
the crossing point and is connected to a glass transition. The 
path for r = 5.50 is close to but below the crossing point and 
close to the A3 -singularity. The third path is connected to a 
gel transition beyond the crossing point. All paths end at an 
A2-singularity given by the respective F. The changes in the 
MSD when approaching the different liquid-glass-transition 
points shall be analyzed using the asymptotic laws for the A 2 - 
singularity in the following. The asymptotic laws for the crit- 
ical relaxation at A2-singularities from Eq. il 11 are compared 
with the full MCT result in Fig.0 For r = 1 .67 (panel A) the 
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FIG. 5: MSD for the SWS at the crossing. Full curves are the re- 
sults for states on the isodiffusivity lines for Dq/D s = 10 5 , 10 7 , 10 10 
marked in Fig.|5| The curves with label c refer to the transition points 
for the value of T indicated. Respective values for the plateaus 6^" 2 
are marked by the symbols +, • and ■ introduced in Fig. [2] In 
the lower two panels, the plateau for the A3-singularity is shown as 
horizontal line. Dotted curves show the leading solution to the crit- 
ical law, (toft)", dashed curves the next-to-leading order for the A 2 - 
singularities, Eq. lilt . Open squares (□) denote the time where the 
solution deviates by 20% from the asymptotic result in Eq. <1 II . An 
effective power law for exponent x = 0.27 appearing at T = 6.63 is 
shown by the dash-dotted line (see text). 

description is similar to that found for the HSS |6]. The expo- 
nent parameter A, A = 0.750 is still close to the one for the HSS, 
X = 0.735. But the time scale t$ = 1.95 differs considerably 
from the value to — 0.425 for the HSS. This is due to a slowing 
down of the dynamics for times where 8r 2 (f) is smaller than 

2 caused by the attractive forces on smaller length scale. The 
exponent for the critical relaxation is a = 0.305. The point 
where the description by Eq. il 11 and the numerical solution 
deviate by 20% of the critical plateau value 6 r c s 2 is marked by 
a square at t w 18 w 9?o. 

Panel B shows the scenario for an approach to an A 2 - 
singularity on the path closer to the A3 -singularity. The ex- 
ponent parameter is increased to X B = 0.857 corresponding 
to a decrease of the critical exponent to a = 0.243. The in- 
creasing importance of the attraction is seen in a decrease of 
the critical localization length representing the plateaus for 
the MSD from 6rf = 0.0318 (labeled by + in panel A) to 
6r c s 2 = 0.0245 (marked by • in panel B). However, the ma- 
jor new phenomenon is the drastic increase of the time scale 
to to ?q = 4 • 10 3 . The critical decay for the A2-singularity 
sets in only for times around t ss 10 6 as indicated by the 
square in Fig. [5] B. There is an additional relaxation pro- 
cess outside the transient ruling the dynamics within the win- 
dow ^ log 10 (f) ^ 4.5. The critical localization length of 
the nearby gel transition yields Sr 2 w 10~ 3 . Therefore, the 
anomalous decay process is not the one related to the gel tran- 
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sition. Rather, it is the decay around the plateau of the close- 
by A3-singularity which appears as a subdiffusive regime with 
almost power-law like variation. This later phenomenon shall 
be explained in detail below. 

In panel C for F — 6.33, the gel plateau is approached with 
t o = 6 • 10~ 3 and the critical relaxation for X c = 0.873 and 
a = 0.232 is described with similar accuracy as discussed in 
panel A. The deviation of 20% is at t = 0.048 = 8 fo and again 
indicated by a square. The comparably large value of X causes 
the leading asymptotic approximation (dotted curve) to devi- 
ate further from the next-to-leading order result. The ampli- 
tude [ZTmsd + K ( a )] m Eq. ( II It is around — 1 in both A and C. 
In this sense, one concludes that the critical dynamics for the 
gel transition is quite similar to the one observed for the glass 
transition. 

The dynamics for the 8r 2 (f) exceeding the respective 
plateaus is quite different for the glass transition shown in 
panel A from the gel transition in panel C. Let us, as usual, 
refer to the process with 8r 2 (f) > 6r£ 2 as an oc-process.The 
a-process shown in panel A is similar to the one in the HSS. 
The crossing of the plateau is followed by a von Schweidler 
relaxation and a crossover to long-time diffusion 1 6 ] . A rescal- 
ing of the time can condense the curves on top of each other, 
a property known as a-scaling. For the dynamics at the gel 
transition shown in panel C, the lower plateau (■) defines 
the onset of the oc-process. The shape of the log8r 2 -versus- 
logf curve differs qualitatively from the one shown in panel A. 
The relaxation around the A3 -singularity plateau causes effec- 
tive power-law behavior with x = 0.27 as shown by the dash- 
dotted line. It is the same phenomenon as observed above 
in panel B. On approaching the gel transition, this subdif- 
fusive regime scales as part of the a-process. This holds if 
the distances to neither the nearby glass transition nor the A3- 
singularity are seriously altered as we further approach the 
gel transition. Under this condition, the A3-singularity and the 
glass-transition singularity influence only the shape of the oc- 
relaxation curves. On the other hand, if the distance between 
the glass-transition and the A3-singularity is changed on the 
path taken, the form of the a-process is also modified. In this 
case, the A3 -singularity is manifested in a violation of the en- 
sealing for the gel transition as found in a recent simulation 
study 1 14]. If the separation from the A3-singularity and the 
glass-transition singularity is sufficiently large, which is true 
for small (p, the dynamics is affected only by the gel plateau 
and directly crosses over from the von Schweidler relaxation 
at the gel plateau to the long-time diffusion. For this reason, 
the exponent x of the effective power law approaches unity 
upon increasing F. 

Figure [6] shows the parameters for the asymptotic descrip- 
tion via Eq. as a function of F along the liquid-glass- 
transition lines for 8 = 0.03. The localization lengths in 
panel A exhibit a jump at the crossing point r° reflecting the 
discontinuous change of f L q . The values for the glass-glass 
transition are also shown down to the A3 -singularity at r°. The 
critical amplitudes /imsd follow the same trend as r L H signaling 
that a change in the localization length also sets the amplitude 
for the relaxation around r£. Panel B shows the two quantities 
in the correction to the critical law. Xmsd shows only small 
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FIG. 6: Parameters for the critical decay at A2-singularities accord- 
ing to Eq. WW . r£ (T) and /jmsd (A) in panel A; k(q) (x) from 
Eq. Kmsd (V), and k(q) + tf M SD (♦) in panel B; and to (0) 
in the panel C. The arrow labeled T° marks the value for the A3- 
singularity, the crossing point. Full and dotted lines are guides to 
the eye to join points on different parts of the glass-transition line for 
< r < and the gel-transition line for T° < T, respectively. 



deviations from the value in the HSS, #msd = — 1 -23. On the 
glass-line at the crossing, .Kmsd = — 1 -57, and on the gel-line it 
reaches Kmsd = — 1.31. At the A3 -singularity, Kmsd = — 1 -64. 
Since away from crossing and higher-order singularities, k(a) 
is always close to zero, the correction to the critical law in 
Eq. il It is dominated by the amplitude A'msd which is neg- 
ative and of order unity there. For this reason, including the 
correction to the critical law in Fig.|5| increases the range of 
applicability considerably in comparison to the leading ap- 
proximation. At higher-order singularities, A, — > 1, and k(a) 
diverges. This is responsible for the increase of the correc- 
tions at the crossing. These corrections change sign when 
K(a) starts to increase. For the case of 8 = 0.03, this happens 
only on the glass-glass-transition line between F° and r°. 

Panel C of Fig.|5]points out the difference in the time scale 
to when coming from small F in the HSS limit or from high 
r, respectively. In the first case, to for the critical law at the 
glass-transition plateau is increasing and eventually diverging 
when the gel transition at the crossing is approached. This is 
because the glassy dynamics of the gel transition determines 
to- For r > r°, fo is orders of magnitude smaller than in the 
HSS since the relevant localization for the gel is encountered 
much earlier in time. On this line of transitions, fo is regular at 
the crossing but diverges at the A3 -singularity. This indicates 
that power laws are an inadequate description of the critical 
relaxation at a higher-order singularity. 

Figure Q displays the parameters quantifying the 
von Schweidler approximation in Eq. H2\ . Panel A 



8 




0123456789 10 



r 

FIG. 7: Parameters for the von Schweidler-law description, Eq. 1 1 21 . 
for 8 = 0.03. Panel A shows the separation parameters a for points on 
the isodiffusivity line for Dq/D s = 10 10 (□ — □). The separation of 
the same points from the A3-singularity, 6i, is shown by the full line. 
The separation £j of points on the liquid-glass-transition for given T 
is shown by filled symbols (♦ • • • ♦), the plus symbol marks 6] for the 
glass-glass transition for T = 5.63. Panel B exhibits the amplitudes 
of the correction in Eq. tl2\ . k(— b) +A^ M sd {<) and k(— b) (x), 
cf. Eq. £5). The values for XjyiSD are me same as shown in Fig. [6] 
Panel C shows the time t~ where the respective critical A2-plateau is 
crossed by the MSD for £>J//> 5 = 10 10 . 



refers to states on the isodiffusivity line D S Q /D S = 10 10 in 
Fig. d The isodiffusivity lines bend away from the crossing 
and this translates into the separation parameters \a\ being 
maximal there. On the same curve, the separation from 
the A3 -singularity | e 1 1 has a minimum around the crossing. 
This also shows that distances in control-parameter space 
as apparent e.g. in Fig. [2] need not necessarily reflect the 
relevant separation parameters of the singularity for the 
asymptotic description. The difference in coordinates of 
the liquid-glass-transition point for F = 5.50 from the A3 
is (A(j),Ar) = (0.085,0.01) while for the crossing point 
(A0,AT) = (0.084,-0.37). This would suggest that the 
former point is closer to the A3 than the crossing point. 
The separation parameters, however, are £1 = —0.028 and 
—0.015, respectively, indicating that the influence of the 
A3 -singularity on the crossing is stronger. Panel B of Fig.0 
displays the correction amplitudes in Eq. d!2i . Kmsd is the 
same as in Fig. [6] and k(— b) shows similar behavior as k(a) 
in Fig. [6] However, as k(— b) is larger than k(o) on the gel 
line it almost compensates the negative values of ZTmsd an( ^ 



^msd + k(— b) is close to zero. 

The time t~ for the onset of the cc-process, i.e. the time 
where the critical plateau is crossed, is shown in panel C. 
When the long-time diffusion is given by the ratio Di/D s = 
10 10 , the plateau in the localization is encountered by the 
MSD for the HSS at t~ = 3 • 10 6 . This is the time when 
the cage around a tagged particle disintegrates and the par- 
ticle starts to diffuse. The increasing attraction for T > in- 
troduces short-ranged bonding among the particles before the 
particles experience the cage. Hence, for the same reason as 
for the increase of fo, this bonding process shifts t~ to higher 
values. When comparing the lower panels of Figs. [6] and 
we observe that for ^ Y ^ 5, the time scales fo and t ~ run 
almost parallel and define a window of six orders of magni- 
tude in time where the cage effect dominates the dynamics. 
For large coupling, T ^ 8, we observe a comparable window 
for the dynamics around the gel plateau, where bonding rules 
the dynamics. Therefore, in both cases the stretching of the 
dynamics is the same what is corroborated by observing that 
X < 0.8 in the mentioned regions l2lll . In this sense, also the 
oc-process of glass- and gel-transition singularities are similar 
if one is unaffected by the other. For 5 < T < 7, or X > 0.8, the 
dynamics is governed by the interference of both mechanisms 
and the emergence of the As-singularity. 
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FIG. 8: Asymptotic description of the MSD near the /^-singularity. 
The full lines are the MSD for states with T = 6.33 and increasing 
(p. Three curves reproduce the results from Fig. [5] C and the last 
one refers to (p = 0.5231. The long horizontal lines show the critical 
plateaus 6rJ 2 for the gel transition at T = 6.33, the /^-singularity 
and the glass transition at the crossing point for T = 5.88. The 
short horizontal lines indicate the corrected plateau 6(r, 2 — P 2 ) for 
the asymptotic laws associated with the respective relaxation. The 
|3-relaxation asymptote around the gel plateau, Eq. , is drawn as 
chain curve labeled p for the solution at D^/D 5 = 10 1() (compare 
text). The chain line labeled vS represents the von Schweidler de- 
scription for the state at (p = 0.5231. For Dq/D s = 10 5 , 10 7 , and 
10 10 dotted and dashed lines show the leading and next-to-leading 
approximation near the A3-plateau in Eq. i 101 . respectively. The 
straight full line labeled x shows the approximation by Eq. I15ai . 
x 7 the corrected power law I16ai . and the dashed line labeled bi the 
approximation by Eq. I16bl . The straight dash-dotted lines show the 
asymptotic long-time diffusion D s t for the respective curves. 
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Figure [8] shows the asymptotic approximation of the oc- 
process for states with F = 6.67 and increasing cp, cf. Fig.|2 
Three plateaus organize the relaxation. First, the gel plateau 
is encountered. Shown here as dash-dotted curve labeled P, 
is the first order description by the full p-correlation function 
from Eq. 0. It continues the description by the critical law 
discussed in Fig. [5] The correction in Eq. (I12> for that A2- 
singularity is close to zero as for almost all gel transitions 
for 8 = 0.03, cf. Fig. This explains why the first-order 
description is so successful in the regime after crossing the 
plateau. After the plateau, the curve for the (3-correlator can- 
not be discerned from the full solution. It extends, acciden- 
tally, also beyond the region of applicability which is limited 
by the A3-plateau. To demonstrate that upon closer approach- 
ing the A2-singularity for the gel transition, the cc-scaling pic- 
ture from Fig. |3 A reemerges, we show an additional relax- 
ation for cp = 0.5231. This has a similar separation parameter, 
O = -10~ 4 , as the curve D S JD S = 10 10 in Fig.|5]A. This last 
curve in Fig.|8]clearly displays the two-step relaxation and is 
described well by the von Schweidler law d!2i . 

The second plateau is associated with the logarithmic relax- 
ation laws. The curvature of the log8r 2 -versus-logf curve is 
positive around the plateau and therefore the leading approxi- 
mation, Eq. d 1 31 , which implies negative curvature, disagrees 
qualitatively. Including the corrections in Eq. dlOt with H(t) 
given by Eq. dl4l >. one gets the dashed lines. These describe 
two decades in time for all curves shown when requiring 5% 
accuracy. The asymptotic laws for the A3 -singularity describe 
approximately half of the relaxation between the gel and the 
glass plateau. In particular, the onset of the effective power 
law discussed in Fig.|5]is captured by the asymptotic approx- 
imation. However, the range of applicability for the logarith- 
mic laws is bound by the neighboring plateaus for gel and 
glass transition. For this reason, the approximations for the 
A3 -singularity do not extend beyond the range shown in the 
figure. In particular, the effective power law with exponent x 
is explained only in the first part by the logarithmic laws and 
is continued by a crossover to the dynamics at the plateau of 
the glass transition. 

To differentiate the effective power law from the power 
laws discussed for the MSD in Ref. l32ll . we show the lat- 
ter for comparison as dotted line in Fig. [8] Let us note first 
that for all states considered we find bj_ > 0. The approxi- 
mation by the leading order power law d!5al > describes one 
and a half decades on the 5% -level as seen for the curve 
Dq/D s = 10 10 . The exponents capture the diminishing slope 
upon approaching the A3-singularity by decreasing from left 
to right, x = 0.331, 0.243, 0.181, 0.163. The corrected power 
law, Eq. ( !16aK yields an exponents' = 0.178 for the last relax- 
ation. This correction comes closer to the effective exponent 
x = 0.27, but improves the description of the effective power 
law only little, as can be seen in the straight full line with 
label x'. When including the curvature bi = 0.0132 in the ap- 
proximation, cf. Eq. J16bi . we find the dashed curve bi, that 
describes the relaxation over three decades in time around the 
A3 plateau. But again it covers only the onset of the effective 
power law. In that sense the effective power law is the ana- 
log of the effective logarithmic decay discussed in connection 



with Fig. 9 of Ref. fl2il . where a crossover from A3- to A2- 
dynamics could explain the observed decay. 

For Dq/D s ^ 10 7 we observe that the curves in Fig.|8]can 
be condensed onto a master curve after the gel plateau. This 
holds for the solutions as well as for the asymptotic approx- 
imations since the distance to the A3 -singularity is no longer 
changed significantly. The decay around the A3-plateau is part 
of the a-process for the gel transition. This a-process con- 
tains also the relaxation around the third plateau in Fig.[8]that 
represents the glass transition at the crossing point. Since the 
distance to this point is relatively large, the asymptotic laws 
are modified by rather large corrections as indicated by the 
plateau correction for the curve labeled D^/D" = 10 10 . De- 
spite the larger distance of the connected glass-transition sin- 
gularity, the last relaxation still slows down the dynamics by 
one decade before the final crossover to the long-time diffu- 
sion. 
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FIG. 9: Variation with the well width. Panel A shows the localization 
length r c s at the crossing point for the glass ( x ) and the gel state (o) as 
a function of 8* — 8 together with the value at the /^-singularity (♦), 
8* = 0.04381. The value of r° s in the HSS is indicated by the arrow. 
Panel B displays the separation parameter — £j (o) and the quadratic 
corrections to the logarithmic relaxation at the crossing point. For 
8 = 8*, 0.03, 0.02, the minimal | e 1 1 is displayed (■) which can be 
reached on the isodiffusivity line D S Q /D S = 10 10 . Panel C displays 
the exponents x (♦), cf. Eq. <15bk and the fit x = 3.05 (8* — 8) as 
dotted line. 

To demonstrate how the crossing scenario in Fig.[8]changes 
when 8 is varied, Fig. |9] exhibits the parameters relevant for 
the description of the relaxation. The three plateaus in Fig. [8] 
are defined by the localization lengths r£. Panel A in Fig. [9] 
shows the variation of the localization lengths. At the A4- 
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singularity, 8 = 8*, all three plateaus join in a single localiza- 
tion length. For 8 < 8*, the localization of a glass state at the 
crossing is larger than the localization of the gel state. This 
difference is becoming more pronounced as 8 decreases. For 
the gel the localization follows 8 and for the glass the localiza- 
tion approaches the value for the HSS. In between there is the 
plateau for the A3 -singularity, which closely follows the local- 
ization for the gel. This limits the regime for the von Schwei- 
dler relaxation after the gel plateau, as observed in connec- 
tion with Fig. [SJ if the A3 -singularity is close. Sufficiently 
far from higher-order singularities, the amplitude in 8r 2 de- 
limited by the localization lengths of gel and glass transition 
exhibits the dynamics defined by a crossover of two different 
A2-singularities. If the A3-singularity is close-by as discussed 
in Fig. [8] logarithmic laws influence the relaxation. 

The influence of the A3-singularity is quantified by the sep- 
aration parameter at the crossing, £j ross , shown for the various 
crossing points in panel B. For smaller 8, the separation in- 
creases and limits the A3-dynamics visible in the relaxation at 
the crossing. The quadratic correction as dominant deviation 
from the logarithmic decay laws is governed by the variation 
of e^ ross while the variation in Kmsd is only small as noted 
earlier [32]. If in an experiment one is limited to a dynamical 
window given by a diffusivity of, say, Dq/D s = 10 10 , this im- 
plies further restrictions to the detection of the higher-order 
singularities. The minimal separation on the isodiffusivity 
curve Dq/D s = 10 10 is shown as in panel B. The expo- 
nent x, cf. Eq. dl5b> . assumed at the crossing point can be 
used as an estimate for the separation from the A3-singularity. 
Since the distance between crossing point and endpoint varies 
quadratically in 8* — 8, cf. inset of Fig. [3] the exponent x at 
the crossing is linear in 8* — 8, cf. Eq. ( I13> . This is shown 
in panel C of Fig. [9] where the exponents can be fitted by a 
linear function. When restricted to the isodiffusivity curve 
Dq/D" = 10 10 , the exponents are larger, accordingly. For 
8 = 0.02 we find x = 0.169 and for 8 = 0.03 the minimal ex- 
ponent is x = 0.095. 



VI. RESULT FOR THE CORRELATION FUNCTION 

The preceding section showed that the dynamical laws at a 
crossing of liquid-glass transition lines can be quite intrigu- 
ing since upon variation of control parameters the separation 
to three different singularities is changed. For the discussion 
of the density correlators <j>o(Y), there enters the wave num- 
ber as a further parameter. Allowing also for a variation in q, 
combines the subtle ^-variation for the logarithmic decay, cf. 
ll32ll . with the g-dependence of the decay at A2-singularities. 
We shall select only a special case which was considered in 
1I2TI1 and found in an experiment lf]~3l fTvIl and also in MD sim- 
ulation m . 

Figure ITOl shows how the dynamics for the states specified 
in the inset is described by the asymptotic laws for differ- 
ent singularities. The interesting feature is the straight line 
piece describing the decay for 0.8 > ty q (t) > 0.6 for state 1 
and 2. This reflects the logarithmic decay caused by the A3- 
singularity. The appropriate plateau value connected with the 




FIG. 10: Logarithmic decay of the density correlation function for 
q = 4.2 near the crossing point for 8 = 0.03. The inset shows part 
of the glass-transition diagram for 8 = 0.03 including the line £1 = 
(dashed). The full curves in the main panel display the solutions for 
states n = 1, 2, 3: (r,cp) = (0.53,5.33), (5.33,0.5361), and (0.53,6) 
which are marked in the inset. Three relevant plateaus are indicated 
by horizontal lines for the gel transition (dashed) at (0.530,6.1) la- 
beled /5 e j, for the A3-singularity (full line) labeled f°, and for the 
glass transition at (0.536,5.33) (short full line) labeled /Sj, ss . The 
plateau values are = 0.954, f = 0.899, and /^ lass = 0.503. Short 

lines show the corrected A3-plateau values f° + 8/ for the three states 
specified. Broken curves show the next-to-leading approximation for 
the logarithmic decay, dotted and dash-dotted curves the leading and 
next-to-leading approximation for the critical decay jSJ in curve 2 at 

fC 

•> glass ' 



A3-singularity is f° = 0.899, and close to the plateau for the 
gel transition f^ el . That the plateaus for gel transitions and 
for the A3 -singularity are close for any wave vector is also 
reflected in the localization lengths in Fig. [9] (A). Therefore 
the logarithmic laws for the A3-singularity have an asym- 
metric range of applicability. The range is rather small for 
shorter times since the gel transition interferes, and consid- 
erably larger for longer times as the critical decay due to the 
glass transition has a more distant plateau. 

The evolution of the dynamics when moving from state 1 
to state 2 is the analog of the dynamics seen in the MSD in 
Fig. |3(B). Only a minor part of the slowing down takes place 
at the gel plateau, the major part from t as 8 to t as 10 4 is 
described by the logarithmic laws around f°. For the solu- 
tions 1 and 2, the approximation by the next-to-leading or- 
der is valid from t as 10 to f « 10 3 and 10 4 , respectively. At 
the A2-singularity for the glass transition, the critical law (jSJl 
is observed. The exponent parameter X = 0.847 implies an 
exponent a = 0.250. The leading f ~ a -law (dotted) describes 
curve 2 successfully for t > 10 6 and adding the correction 
(dash-dotted) improves that range by almost two decades. 
Curve 2 demonstrates how different asymptotic expansions 
complement one another: The logarithmic laws describe the 
decay from above f° down to §q(t) > 0.7 and Eq. (jSJi approx- 
imates successfully the region from tyq(t) < 0.7 to the critical 
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plateau /f lass . That the slope of the decay becomes smaller 
below f° is a clear indication of a closer approach to a higher- 
order singularity, as prefactor B in Eq. Jl 31 vanishes with the 
square-root of the distance from the A3-singularity. 

When taking another path from 1 to state 3, the distance to 
the A3-singularity remains largely unaltered and we find the 
counterpart of Fig.[8]for the MSD. The dynamics is ruled by 
an approach to the gel transition and the complete decay be- 
low /g el is part of the a-process. This a-process for the gel 
transition scales by a shift along the respective plateau /| el 
with only minor deviations due to changing separations to the 
glass-transition line and the A3 -singularity. No clear two-step 
process is observed for curve 3 for two reasons. First, the A2- 
dynamics below f^ ei is limited by the logarithmic laws for the 
A3-singularity. Second, the complete decay seen in curve 3 
requires more than ten decades, but only for t < 10 2 the decay 
takes place above the plateau /| el . Hence, the decay onto the 
plateau is too close to the transient dynamics to exhibit a clear 
critical decay. Moreover, the exponent parameter in the vicin- 
ity of the A3 -singularity is already rather high, X = 0.89, so 
the critical law t~ a is stretched considerably. As in the MSD 
shown in Fig. [8] for the last curve, moving closer to the gel 
transition, the two-step process typical for an A2 -singularity 
reemerges. 

VII. CONCLUSION 

The relaxation scenarios for line-crossings near higher- 
order glass-transition singularities were presented in this 
work. Three different singularities influence the dynamics in 
that region of the glass-transition diagram, and asymptotic ex- 
pansions around each of these are necessary to successfully 
describe the complete relaxation patterns. Each singularity 
is associated with a characteristic plateau value as shown for 
the localization lengths for the MSD in Fig. [9] A and for the 
Debye-Waller factors fS for $q(t) in Fig.^| The position of 
the different plateau values arranges the successive steps for 
the relaxation in time. 

The plateau of the gel transition is encountered first. It is 
approached by the relaxation with the critical law of the A2- 
singularity, cf. Fig.[5]C. The dynamics after crossing the gel 
plateau is described by the von Schweidler law related to the 
A2-singularity for the gel transition, cf. Fig- (HI before the log- 
arithmic laws at the A3 -singularity become valid. The latter 
have been studied extensively and imply a subdiffusive power 
law with exponent x, cf. Eq. dl5b> for the MSD at specific 
points in control-parameter space where bi in Eq. J16bl > van- 
ishes [:32f] . However, for a region near the crossing where 
£>2 > 0, an effective power law with exponent x can be identi- 
fied in Fig.|5]C. The onset of this behavior is described by the 
asymptotic dynamics around the A3 -singularity while the ex- 
tension to later times originates from a crossover to the critical 
dynamics at the plateau of the glass transition, cf. Fig. [8] Both 
the asymptotic power law fl^l and the crossover scenario flUl 
have been found for the MSD in recent computer simulation 
studies. 

A similar crossover which yields the f x -relaxation in the 



MSD is responsible for an effective logarithmic decay in the 
correlation functions for wave vectors that are accessible in 
typical light-scattering experiments. Again, the dynamics be- 
tween the plateau for the A3-singularity and the plateau for the 
glass transition assumes a variation linear in lnf, cf. Fig. 1101 
Most of this behavior is fitted satisfactorily by two different 
asymptotic laws and is therefore clearly differentiated from 
the asymptotic logarithmic decay at higher-order singularities 
which is expected only for large values of the wave vector 
l32ll . Nevertheless, also the effective logarithmic decay can 
serve as a clear signature of a line crossing and hence for the 
existence of higher-order singularities. The decay analyzed 
in Fig.[H)|has been identified as a typical scenario in systems 
with short-ranged attraction in experiment fl3ll . theory [21], 
and computer simulation 1131 . 

The last relaxation step of the complete decay in the vicinity 
of the line crossing occurs at the plateau for the glass transi- 
tion and is similar to the scenario known from the HSS as seen 
in Fig.|5]A. Only after having crossed this last plateau, the dy- 
namics enters the long-time diffusion limit. Each of the relax- 
ation steps discussed above can be more or less pronounced 
depending on the separation from the related singularity in 
control-parameter space. It can be inferred from Figs. |5J 
and [9] that in a certain region around the higher-order singu- 
larities, the presence of the latter singularities introduces large 
corrections to the asymptotic laws at gel- and glass-transition 
points. Outside this region, however, the use of the conven- 
tional A2-scenario is justified and the asymptotic approxima- 
tion varies only little there. Hence, the dynamics near any 
state on the entire surface of liquid-glass and liquid-gel transi- 
tions can be characterized by the parameters of the asymptotic 
approximations. 

The variation of the final long-time diffusion can be used 
to map the theoretical glass-transition diagram to the experi- 
mental control-parameter space and thus locate higher-order 
glass-transition singularities at least approximately. The map- 
ping proposed in this work could be used to estimate the lo- 
cation of an A3-singularity in 1 15] by extrapolation, cf. Fig. 
|4] and facilitated the identification of an A4-singularity in a 
recent computer simulation study |16]. Within the Percus- 
Yevick approximation for S c/ , the A3 -singularities are behav- 
ing similar to the critical points of the fcc-fcc binodal 1 40] : 
Upon changing the well width 8, MCT endpoints and critical 
points vary only little in the attraction strength F as seen in 
Fig. [3] When using the structure factor in mean-spherical ap- 
proximation, this behavior is different. But this difference is 
eliminated after identifying the glass-transition diagrams for 
both closure relations at the crossing points. For 8 = 0.03, the 
densities of endpoint and critical point are in accord reason- 
ably, while the higher value for T fixes the A3 -singularity in 
the metastable region with respect to the isostructural phase 
transition. 
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APPENDIX A: CONSISTENCY OF THE NUMERICAL 
SOLUTION 

Glass-glass-transition points and higher-order singularities 
were calculated for Figs. ^ EE] and |9] The expeditious and 
accurate identification of these singularities is also crucial for 
the evaluation of the asymptotic approximations. Therefore, 
some notes concerning the numerical solution of Eq. 10 shall 
be discussed in this appendix. For the determination of liquid- 
glass-transition points a robust method of nested intervals can 
be applied anticipating the jump from zero to a finite value 
in the glass-form factors f„ at the respective A2-singularity. 
This procedure works also at an A4-singularity which is also a 
liquid-glass-transition point. For a glass-glass-transition point 
the discontinuity in the glass-form factors takes place between 
finite values and the jump in the f q becomes smaller when ap- 
proaching the A3 -singularity and observing a discontinuity in 
the glass-form factors becomes increasingly difficult. There- 
fore a different criterion shall be used. To this end, coeffi- 
cients from the expansion of the RHS of Eq. (|2ji are required, 

cf. E3. 

4lU = ^-ftirr^ffiflfc-afO (A1) 

x(l-£)- ■■(!-.£). 

At a glass-transition singularity, Eq. is no longer invertible 
which is signalled by the maximum eigenvalue E of the so- 
called stability matrix A^ c approaching unity from below @|. 
The evolution of E in the vicinity of an A2-singularity is given 
by a square-root in some control parameter v, 1 — E « \/v — v c , 
for the strong-coupling side v > v c . Monitoring the eigenval- 
ues can be done with high precision and allows for an extrap- 
olation in control parameters which can reduce the numerical 
effort considerably. At an A3 -singularity, the eigenvalue is ap- 
proaching unity from either side on generic paths in control- 
parameter space through the singularity. The variation is given 
by 1 — E °< (v — v°) 2 / 3 which follows from generic properties 
of the singularity 14411 . 

It is clearly seen in Fig.[^that at a glass-glass transition 
only the eigenvalues for the strong coupling side, (p > (p e , go 
to unity and follow the square-root law. At a liquid-glass tran- 
sition the eigenvalues for <p < <p c would be zero, however, in 
the glass due to continuity they are finite, smaller than unity 
and jump to a critical value only at the glass-glass-transition 
points. For the A3-singularity this discontinuity vanishes and 
the eigenvalues show the variation with the power 2 /3 on both 
sides of (p°. The deviation from that law for larger distances 
with cp < cp° is due to the increase of the eigenvalues at the 
liquid-glass transition at <p = 0.540693. Deviations close to 
the A3-singularity on the other hand indicate the precision 
of five digits in the control parameter <p for the determina- 
tion of V°. The deviation of E c from unity is a measure for 
the accuracy of the critical points. In this work a value of 
1 — E c < 10~ 3 was assured for all the transition points shown 
in this work. 




log,o lt P- ( P C| 



FIG. 11: Eigenvalues E upon approaching a glass-glass transition 
for 8 = 0.02, r = 7.75, and (p e = 0.540965015. The deviation from 
unity, 1 — E, is shown for (p < (p f (open circles) and for (p > (p f (filled 
circles) together with the square-root ^/|(p — (p c | (dashed). The corre- 
sponding eigenvalues for the A3-singularity at 8 = 0.02, T° = 6.646, 
and (p° = 0.5680321 are denoted by open squares for (p < (p° and 
by the filled squares for (p > (p° . The full line shows the power law 
|cp-<p°|2/3. 
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FIG. 12: Parameter fi 2 = 1 - A, in the SWS for 8 = 8* and 0.03. 
Values for liquid-glass transitions are shown as full lines, for glass- 
glass transitions as filled circles. The dashed lines show the laws 
(r-T*) 2 / 3 for the A4- singularity and /J2 x (r-T ) 2 / 3 for the 
A3-singularity. The squares indicate a deviation between result and 
approximation of 5%. 



Despite being useful as an extrapolation scheme, the 
generic laws close to the singularities can also serve as consis- 
tency check for the numerical results. This was already shown 
in the inset of Fig.|5]for the distance of the crossing point from 
the A3-singularity. There, the control parameters close to the 
A4-singularity were related in a quadratic polynomial. As an- 
other quantity we utilize the exponent parameter X which ap- 
proaches unity at higher-order singularities. /J2 = 1 — A. is also 

given by coefficients from Eq. iAl\ . ^2 — 1 — a *q^qk^k 2 ak i ak ^^ 
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where summation over repeated indices is assumed and a* and 
a q denote the left and right eigenvectors of the stability matrix 

A^ c , respectively. 

Figure shows that close to higher-order glass-transition 
singularities the exponent parameters X = 1 — /i2 calculated 
numerically obey the asymptotic approximation by the re- 
spective power laws with reasonable accuracy. For the Aj- 
singularity the description works down to A, = 0.85 and in- 



cludes both glass-glass transitions and liquid-gel transitions. 
The Azt-singularity is described by the asymptotic law for 
X ^ 0.93 on the line of gel transitions and for X ^ 0.9 on 
the line of glass transitions. The exponent parameters for dif- 
ferent potentials fall on top of each other close to their A4- 
singularities jU That the asymptotic approximation is ap- 
plicable for a similar range in control parameters underlines 
the universality of the A4-singularity. 
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